Hyperfine meson splittings: chiral symmetry versus transverse gluon exchange 
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Meson spin splittings are examined within an effective Coulomb gauge QCD Hamiltonian incorpo- 
rating chiral symmetry and a transverse hyperfine interaction necessary for heavy quarks. For light 
and heavy quarkonium systems the pseudoscalar-vector meson spectrum is generated by approxi- 
mate BCS-RPA diagonalizations. This relativistic formulation includes both 5* and D waves for the 
vector mesons which generates a set of coupled integral equations. A smooth transition from the 
heavy to the light quark regime is found with chiral symmetry dominating the -K-p mass difference. 
A reasonable description of the observed meson spin splittings and chiral quantities, such as the 
quark condensate and the tt mass, is obtained. Similar comparisons with TDA diagonalizations, 
which violate chiral symmetry, are deficient for light pseudoscalar mesons indicating the need to 
simultaneously include both chiral symmetry and a hyperfine interaction. The rjb mass is predicted 
to be around 9400 MeV consistent with other theoretical expectations and above the unconfirmed 
9300 MeV candidate. Finally, for comparison with lattice results, the J reliability parameter is also 
evaluated. 
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I. INTRODUCTION 

The hyperfine interaction has a long and distinguished 
history beginning with the hydrogen atom where it cor- 
rectly describes the transition responsible for the famous 
"21-centimeter line" in microwave astronomy. Taking the 
non-relativistic reduction of the one-photon-exchange in- 
teraction, the hyperfine potential has the form 



Vhvv — const — 2- • — '^5(r) 



(1) 



for particles of mass Mi and spin cr^. This potential 
gives an accurate description of the triplet-singlet split- 
ting in positronium. When implemented in the simple 
additive quark model having meson mass M and con- 
stituent quark masses M.q=u.d,s , 



M ^Mn+Mn 



q Mq 



(2) 



it produces a remarkably good description of the light 
meson spectrum using A k, 160 Ai'^ MeV. However, to re- 
produce the splitting in charmonium, using M.q=c ~ 1500 
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MeV, requires a hyperfine strength of at least AA, while 
a similar model for baryons uses a much weaker value, 
approximately A/ 3 P|. Hence attempts to comprehen- 
sively describe spin splittings in hadrons with a simple 
hyperfine interaction leads to over an order of magni- 
tude variation in the potential strength. Although this 
may merely reflect the simplistic nature of the additive 
quark model, more extensive models also have difficulty 
in obtaining a consistent description for both mesons 
and baryons with the same hyperfine interaction. Pos- 
sibly related, other problems arise when considering the 
hyperfine interaction in quark model applications. For 
example, the non-relativistic reduction of the one gluon 
exchange interaction between quark pairs gives rise to hy- 
perfine, spin-orbit and tensor interactions at order 1/A^^. 
Unfortunately this structure does not describe heavy me- 
son spin splittings well and must be supplemented with a 
spin-orbit term which is argued [3] to emerge from the re- 
duction of a scalar confinement potential. Although this 
prescription is reasonably successful for heavy mesons, 
it generates too large spin-orbit splittings in baryons Q 
and the sign of the scalar confinement potential must be 
inverted 4] . 

Further clouding this matter is the role of chiral sym- 
metry. Because the pion is regarded as the Goldstone 
boson of broken chiral symmetry, the ir-p mass differ- 
ence should be dominated by the same nonperturbative 
dynamics producing spontaneous chiral symmetry break- 
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mg |S|,|y,|Y|,|^. This observation is in conflict with the hy- 
perfine splitting in typical constituent (non-chiral) quark 
models which apply Eq. to light hadrons. Indeed 
in such models the hyperfine potential plays a dual role 
of generating spin splittings and producing a very light, 
"chiral" pseudoscalar meson. Because of the quite large, 
over 600 MeV, ir-p and 400 MeV K-K* mass differences, 
an appreciable hyperfine splitting is required. This in 
turn makes it necessary to use a more complex interac- 
tion, rather than a simple potential, to simultaneously 
describe hadrons not governed by chiral symmetry, such 
as excited state light mesons, heavy mesons and baryons 
which all have smaller spin splittings, roughly from 50 
to 300 MeV. As discussed in Ref. the nonperturba- 
tive aspects of spin splittings in light quark hadrons, es- 
pecially those involving pseudoscalar mesons, clearly re- 
quires a more sophisticated model treatment. 

The purpose of this paper is to address the above 
issues and to provide a deeper understanding of me- 
son spin splittings by examining the hyperfine interac- 
tion in a theoretical framework which incorporates chi- 
ral symmetry. A related goal is to also provide an im- 
proved hadron approach embodying many of the fea- 
tures of QCD with a minimal number of parameters 
(i.e. current quark masses and one or two dynamical 
constants). This formulation is based on the Coulomb 
gauge Hamiltonian of QCD and approximate diagonaliza- 
tions using the Bardeen-Cooper-Schrieffer (BCS), Tamm- 
Dancoff (TDA) and random phase approximation (RPA) 
many-body techniques. These methods have been pre- 
viously applied to chiral symmetry breaking [H|. glue - 
balls |ll|l5,'l3i, hybrids and mesons [SEElIi- 

The renormalization group methodology has also been 
applied0i |20| to improve this formalism. Finally, this 
work incorporates an effective QCD longitudinal confin- 
ing potential Q along with a generalized version of the 
transverse hyperfine interaction employed in an earlier 
hyperfine study |2ll |. 

An important aspect of the following discussion is that 
the random phase approximation is capable of describing 
chiral symmetry breaking. Indeed, the RPA pion mass, 
Mtt, satisfies the Gell-Mann-Oakes-Renner relation dic- 
tated by chiral symmetry 



2m„ 



m) 



P 

J TT 



1/2 



(3) 



where ruq is the current quark mass, is the pion decay 
constant and {qq) = (0|^^'|0) is the quark condensate. 
In contrast, the TDA does not respect chiral symmetry 
(the RPA, but not TDA, meson field operator commutes 
with the chiral charge 0, ^3)- Comparing RPA and 
TDA masses therefore permits a quantitative assessment 
of the relative importance of both chiral symmetry and 
the hyperfine interaction for the 7r-p mass splitting. A 
key result of this work is that chiral symmetry domi- 
nates this splitting and accounts for about 400 MeV of 
the mass difference. Another is that with the same inter- 
action, the mass differences between the radially excited 



TT and p states are also reproduced, as well as the pseu- 
doscalar and vector states in charmonium and bottomo- 
nium, thereby demonstrating the universality of this ap- 
proach. Lastly, the inclusion of the hyperfine interaction 
improves the model description of the quark condensate. 

In the next section the model Hamiltonian is specified 
and the BCS and RPA equations are formulated. Section 
III presents numerical results and details sensitivity to 
different hyperfine interactions. Finally, conclusions are 
presented in section IV. 



II. EFFECTIVE HAMILTONIAN AND 
EQUATIONS OF MOTION 

A. Model Hamiltonian 

As discussed above, the model Hamiltonian is taken to 
be that of Coulomb gauge QCD. The Coulomb potential 
is evaluated self-consistently in the mean-field gaussian 
variational ansatz 4]. The general form of this effective 
Hamiltonian in the combined quark and glue sectors is 

H„ff = H, + H, + H,g + Vc (4) 

Hq = j dx«'t(x)(-iQ; • V -t- TO,/3)^'(x) (5) 

Hg = Tr j dx[n"(x) • n°(x) + B°(x) • B°(x)](6) 

Hqg = g ydxJ'^(x).A"(x) (7) 

= -\J dxdyp'^(x)l/(|x-y|)p'^(y) . (8) 

Here g is the QCD coupling, 5* is the quark field, A° are 
the gluon fields satisfying the transverse gauge condition, 
V • A° = 0, a = 1, 2, ...8, n° are the conjugate fields and 
B'' are the nonabelian magnetic fields 



B" ^ V X A" + ig/° 



"A" X A'^ 



(9) 



The color densities, p°(x), and quark color currents, J", 
are related to the fields by 

p"(x) = «'t(x)r°«'(x)-i-r''"A''(x)-n"(x) (10) 

J" = ¥{x)aT''^{x) , (11) 

where T" = ^ and /"^"^ are the SU3 color matrices and 
structure constants, respectively. 

Since this work focuses on the quark sector Hg is omit- 
ted and the quark-glue interaction, -ffgg, is replaced by an 
effective transverse hyperfine potential, Vr, discussed be- 
low. The choice for the longitudinal Coulomb potential, 
Vc, then completes specification of the model. 

To lowest order in g, the Coulomb gauge potential V 
is simply proportional to 1/r. As is well known, in a 
few body truncation this is insufficient for confinement 
and fails to reproduce the phenonienological meson spec- 
trum. Using the Cornell potential for V resolves these 
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issues but produces an ultraviolet behavior necessitat- 
ing the introduction of a model momentum cutoff. In- 
stead, an improved dynamical treatment Q is adopted in 
which both the gluonic quasiparticle basis and the con- 
fining interaction were determined self-consistently and, 
through renormalization, accurately reproduced the lat- 
tice Wilson loop potential. The resulting interaction has 
a renormalization improved short ranged behavior and 
long-ranged confinement. It is similar to the Cornell po- 
tential and has a numerical representation in momentum 
space that is accurately fit by the analytic form 



C{P) 



1.07 



log 



-0.62 



hO.82 



r + 1.41 



L{p) 



12.25 m„ 



for p > nig 

for p < nig . 

(12) 

The low momentum component is numerically close to 
a pure linear potential, L ~ —Sira/p'^. The other term 
represents a renormalized high energy Coulomb tail. The 
only free parameter is nig « 600 MeV which sets the scale 
of the theory and is equivalent to a string tension. 

Both the exact and model QCD Coulomb gauge Hamil- 
tonians do not explicitly contain a hyperfine-type inter- 
action, however, perturbatively integrating out gluonic 
degrees of freedom generates a quark hyperfine interac- 
tion with the form a.i-a2- One can also formally generate 
this Lorentz structure using Maxwell's equations to sub- 
stitute for the gluon fields in Eq. (7). More generally, 
in the Hamiltonian formalism an effective hyperfine in- 
teraction arises from nonperturbative mixing of gluonic 
excitations (such as hybrids) with the quark Fock space 
components of a hadron's wavefunction. However, the 
cxi ■ CX2 Lorentz structure is expected to persist and 
Ref. obtains a 0.1-0.2 hyperfine potential with specific 
spatial form using the linked cluster expansion method 
to eliminate hybrid intermediate states. 

Because contributions from gluonic excitations and hy- 
brid states are difficult to calculate, this work studies sev- 
eral different parameterizations of the following generic 
transverse hyperfine interaction 



Vt 



1 



dxdyJf(x);7.,(x,y)Jf(y) 



where the kernel Ua has the structure 



(13) 



(14) 



reflecting the transverse gauge condition. In specifying 
the potential U it is useful to realize that perturbatively 
U — > as/|x — y| where as — /Att. Also, the form of 
the one gluon exchange potential and the nonperturba- 
tive mixing with hybrids makes it clear that the hyperfine 
kernel should not include a confining term. This is impor- 
tant for infrared divergent gap equations . Consistent 
with these points, the following four kernels are utilized 
and numerically compared to document hyperfine model 
sensitivity. 



The non-relativistic quark model advocates a regulated 
contact interaction Thus model 1 is a simple square 
well interaction defined by 



Ui{p) 



for p > A 

-Uh for p < A 



(15) 



with strength, Ut, and range, A. 

Model 2 is a variation of a pure Coulomb potential 
(reflecting a transverse zero mass gluon exchange) and is 

TT ( \ / forp> nig 

U2{P) = [_QU iovp<nig. (16) 

Similarly, model 3 incorporates a modified Coulomb 
potential corresponding to an ultraviolet Coulomb tail 
matched to a constant in the infrared. This potential is 



_ J C{p) for p> nig 
-Ch for p < m 



9 ■ 



(17) 



Finally, model 4 is a Yukawa-type potential corre- 
sponding to the exchange of a constituent gluon with 
a dynamical mass. This is given by 



j C{p) for p>mg 
Uiip) = <! _^£Z^ forp<mg 



(18) 



For the latter three models, the Coulomb potential 
C{p) is the same as in Eq. ifT^ and the constant, Ch, 
is determined by matching the high and low momentum 
regions at the transition scale nig . In this analysis several 
different matching points (e.g. p = mg,2mg,Smg) for a 
given transition scale mg were numerically examined but 
no qualitative differences were found. 

In the following, angular integrals are denoted by 



Vn{k,q) = J dxV{\k~q\)x" 



(19) 



where a: = k • q. The auxiliary functions 

>F(|k-q|).^(|k-q|) -(^^+f)-^,y + -^) (20) 
and 

1 - 1-2 

(21) 



Z(|k-q|)^C/(|k-q|)^L^ 



|k-qP 

are also introduced which arise from the operator struc- 
ture of Eq. ifTH) . 



B. Gap equation 

Calculations are most conveniently made in momen- 
tum space with constituent quark operators. These are 
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obtained by a Bogoliubov transformation (BCS rotation) from the current quark basis to a quasiparticle quark basis 
represented by particle, B, and antiparticle, D, operators 



(22) 



where ii is a color vector with i — 1, 2, 3 and A denotes helicity. The Dirac spinors are functions of the Bogoliubov 
gap angle 0fe and can be express in terms of the Pauli spinors x\ 



V- 



kA 



1 

1 

71 



VI + sm(t>k XX 
VI - sin^fc cr ■ k xa 

-VI - sin0fc cr ■ k i(T2X\ 



(23) 
(24) 



VI + sincjjk i(72X\ 

Note the additional factor ia2 in the spinor V-kA which differs from the convention used in Refs. 0, 0|. The 
Bogoliubov angle can be related to a running quark mass, Mq{k), and energy, E{k) — ^ M^{k) + k^, by 

M,(fc) 



Sk 



Ck 



sm (pk 



cos (pk 



E{k) 
k 

m 



At high k, Mq(k) nig, while for low k a constituent quark mass can be extracted, A4q = Mq{0). 
Minimizing the vacuum energy with respect to the gap angle yields the mass gap equation 



k Sk - m„ Ck = 



dq 

12^ 



[{skCqX - SqCk)V{\\<. - q|) - 2ckSqU{\k - q|) + 2cqSkW{\k - q| 



which, after angular integration, reduces to 

k Sk - mq Ck ^ 
Finally, the quasiparticle self-energy is 

ek = ruq Sk + k Ck 



oo 2 



q dq 
67r2 



oo 2 



?kCq{Vl + 2Wo) - SqCkiVo + 2Uo)] 



q dq 
6^ 



[skSqiVo + 2Uo) + CkCqiVi + 2Wo)] 



(25) 
(26) 

(27) 
(28) 
(29) 



C. Meson RPA equations 

The pion RPA creator operator, XapB^^Dj^ — yapBaDp, contains two wavefunctions with color a, 6, spin a, /9, 
momentum, k, and radial, v, indices given by 

A'^(k) = -^'S^Ill^^-^X-ik) (30) 
V4^ V2 V3 

3;-(k) = 1 !(^%F-(fc) . 
V4^ 72 V3 



Diagonalizing the effective Hamiltonian in the RPA representation yields two coupled radial equations valid for any 
equal quark mass pseudoscalar meson 

/■oo 2 J 

2ek X^k) + \^ {K{k, q)X'^{q) + K\k, g)y-(g)] = Mi;X''{k) (31) 

2ekY''{k)+ r l^[K{k,q)Y''{q) + K'{k,q)X-^{q)] = -M^Y^k) , 
Jo Stt^ 



with kernels 



K{k,q) = {l + SkSq)Vo + 2{l-SkSq)Uo+CkCq{Vl~2Wo) (32) 
iV(fc, q) = (1- SkSq)Vo + 2(1 + SkSq)Uo - CkCqiVl - 2Wo) . (33) 
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The above equations yield a zero mass pion in the chiral Hmit 
self energy equations to obtain 



|24| as demonstrated by combining the gap and 



Sfeefc = rriq 



q^dq 
6^ 



and then substituting this expression into the RPA equations, first multiplied by Sk- For mq — 0, X(k) and y{k) 
become proportional to s^. This immediately yields the eigenvalue M^r = in accord with Goldstone's theorem. We 
have numerically confirmed this to a precision of 5 keV in solving the RPA equations. 

Constructing the p and other vector meson RPA wavefunctions requires three spin projections. Also, both S and 
D orbital waves now contribute and the general solution is 



(34) 
(35) 



It is reasonable to assume exact isospin symmetry (degenerate quark masses, m„ — rrid), which permits suppression 
of this quantum number. The four wavefunctions having unit norm are then 



4^ \/2 V3 

1 3,.~ .~ 1 Sab 



x:{k) 



(36) 



_-(k.<Tk--^)za2-2=X,^(fc) 



1 (Tia2 Sal 

4^ \/2 



x: = 

Exploiting the symmetry of the RPA kernels, under transposition and simultaneous k q exchange, reduces 
the number of independent kernels to six: Kj^j^, ^xx^ ^'xx' ^xy^ -^xV' ^xy- '^^^ other required kernels 
can be obtained from these using Kyy = Kxx, Kxy = Kyx for all four angular momentum combinations and 
K'^''{k,q) = K'"^{q,k) for all four X-Y combinations. The p RPA equations are 



i=^(k..k-ia)za.^r,''(fc) 















\Yd 



q^dq 
6^ 



T^ss Ti'sd K'ss irsd 

^XX ^XX ^XY ^XY 

J^ds T^dd T^ds T^dd 

^XX ^XX ^XY ^XY 

jy-ss irsd T^ss TTsd 

^YX ^YX ^YY ^YY 

T^ds T^dd r^ds T^dd 

^YX ^YX "-YY "-YY 



Yr 



Ml: 



(37) 



and the integration is performed after multiplication with the column wavefunction vector. 
Finally, the six independent kernels are 



K 



XX 



^XX 



i [3(1 + .Sfe)(l + s,)Vb + (1 - sq){l - Sk){AV2 - Vo) + 2cfcC,(31/i + 2C/i + 2kqZo) + 
2(1 + Sk){l - Sg)(-i7o + 2fc2Zo) + 2(1 - Sfe)(l + Sq){~UQ + 2q^Zo)] 
V2 



K 
K 



XY 

d 
XY 



= ^ [(1 - Sq){l - Sk){V2 - Vo) + (1 + sk){l - s,)(-2;7o + fc'Zo) + 2cuCq{2Ui 

(l_5,)(l + sJ([/„_3[/2+g2^0)] 
= \ [-(1 - SkSq)Vo - (1 + Sk){l + Sq)U^ + (1 - Sk){l - Sq){Uo - 2U2) + CkCq{Vl - 

= ^ [2(1 + 5fc)(l - Sq)Vo + (1 + Sq){l - Sk){3V2 - Vq) + (1 + Sg){l + .Sfc)(-2C/o 



kqZa) + 
- 2Ui + 2kqZn) 



(38) 

(39) 
(40) 

(41) 



J<xx 



K 



dd 
XY 



+ (1 - Sfe)(l - Sq){-Ua - C/2 + 3fc2Zo) - 2cfcC,(2V\ + 2{/i + kqZ^)] 
^ [3(1 + Sfc)(l + Sq){W2 - Ml) + (1 - Sfc)(l - Sq){V2 + 5^)) + 4cfcC,(31/i + + kqZ^) + 
(1 + Sfc)(l - Sq){-Wo - 3U2 + {k^ + 9q^)Zo) + (1 - .sfe)(l + Sq){-5Uo - SU2 + i9k^ + q^)Zo)] (42) 
^ [2(1 ~ SkSq)i3V2 - Vo) - AckCqiVi + 4C/i - kqZo) + (1 + Sfc)(l + s,)(t/o - 9C/2 + 3ik^ + q^)Zo) 
+ (1 - .H.)(l - Sq){-7Uo ~U2 + 3(fc2 + q^)Zo)] . (43) 
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III. NUMERICAL RESULTS AND DISCUSSION 

A. Meson spectra 

The numerical techniques for solving the gap equation 
and dia gona lization for the meson eigenvalues are given 
in Refs. The four hyperfine model inter- 

actions were each determined by fitting the charmonium 
Vc-J/"^ splitting. In addition to adjusting the potential 
parameters, Uh, A for model 1 and nig for models 2, 3 and 
4, it was also necessary to significantly reduce the current 
charm quark mass from values typically used which is dis- 
cussed below. The hyperfine potentials are summarized 
in Table n 

TABLE L Different hyperfine efi'ective interactions. 



Model 


Parameters 


1. square well 


Uh = 5GeV-^ A = 3mg = 1.95GeV 


2. Coulomb 


TUg = 0.43 GeV 


3. modified Coulomb 


mg = 0.6 GeV 


4. Yukawa 


rrig = 0.6 GeV 



Using these potential parameters the remaining light 
and heavy pseudoscalar and vector meson spectra were 
then predicted. The 7r,p, 77^, J/^*, rji, and T ground and 
excited states are listed in Tabled While all four models 
provide similar, reasonable meson descriptions, potential 
4 emerges as the preferred model. Note that it was again 
necessary to reduce the current quark masses. 

It is significant that the RPA Hamiltonian approach, 
using any of the four hyperfine interactions, can simulta- 
neously describe both the large ir-p mass difference and 
the small ir'-p' and charmonium splittings. Figure ^ fur- 
ther illustrates this by comparing the RPA (solid circles), 
TDA (squares) and observed (diamonds) hyperfine split- 
tings versus the spin averaged pseudoscalar and vector 
meson mass. Notice that similar to observation both the 
RPA and TDA yield a rapid decline in the spin split- 
tings with increasing meson mass, however only the RPA 
can describe the sizeable ir-p difference. This is because 
the RPA consistently implements chiral symmetry while 
the non-chiral TDA predicts a pion mass that is too large 
(about 500 MeV). Similar to the findings of Ref. [ll|, chi- 
ral symmetry is clearly the dominant effect in the large 
TT-p splitting, accounting for almost 70% (roughly 400 
MeV) of the mass difference. 

It is insightful to contrast this with non-relativistic 
quark model treatments which use a constituent quark 
mass about half the multiplet average and describe the 
splitting as a l/Ai"^ dependence characteristic of rela- 
tivistic corrections. Because the quarkonium Bohr or- 
bit scales inversely with quark mass, it is possible for 
these models to describe both light and heavy meson 
splittings with the same short ranged hyperfine poten- 
tial although it does require tuning of parameters. In- 
deed as Ref. [21 details, a comprehensive meson descrip- 



TABLE II: Calculated masses, condensates and data in MeV 
(rounded to the nearest 5 MeV). 



Quantity 



''rrib 
Mu = Md 

Ma 

Mb 

- < qq>3 



7r(1300) 
(1800) 



M, 
M, 



Mp 

p(1450) 



p{170D) 



M, 



M 



7)c(is) 

rjc(2S) 
rjc(3S) 



J/*(1S) 



*(2S) 



*{3770) 

Mr 

Mt(2s) 



Model 1 Model 2 Model 3 Model 4 experiment 



1 


1 


1 


1 


1.5 - 8.5 






oou 


"ii n 

U lU 


1 nnn 1 /ton 


OOOU 


ZoUU 


Z i OU 


971 n 

Z / lU 


4UUU-40UU 


00 


140 


00 


1 nn 
lUU 


zUU - oUU 


luyu 


lloU 


luyu 


1 non 
luyu 


lOUU 


4Uzo 


/inon 
4UzU 


oyoU 


oyDO 


40UU-01UU 


150 


200 


165 


180 


220-260 


195 


150 


190 


190 


138 


1430 


1150 


1350 


1370 


1300 


2170 


1650 


2085 


2100 


1801 


820 


755 


780 


795 


771 


1480 


1150 


1405 


1420 


1465 


1725 


1305 


1605 


1620 


1700 


2980 


2950 


2990 


2985 


2980 


3660 


3400 


3615 


3625 


"3631 


4210 


3720 


4090 


4100 


7 


3110 


3130 


3110 


3130 


3097 


3740 


3470 


3670 


3680 


3686 


3780 


3490 


3685 


3695 


3770 


9395 


9360 


9415 


9395 


? 


9465 


9440 


9470 


9460 


9460 


9915 


9705 


9880 


9870 


10023 



^adjusted 

"The Belle collaboration ^ reports 3654 ± 0.014 

Hyperfine splitting vs. average multiplet mass 



• RPA 

B TDA 

<i> EXP/NRPQCD 



i*. 

urn* 



4000 6000 
(3M(l--)+M(0-+))/4 



•f 

I ,ni B 

8000 10000 



FIG. 1: Hyperfine splitting versus spin averaged meson mul- 
tiplet mass. 



tion can be obtained by using a potential with a com- 
plicated, mass dependent short ranged smearing. Al- 
ternatively, the RPA-BCS formalism, which dynamically 
generates a constituent mass, reproduces this behavior 
via chiral symmetry and a simpler, weaker hyperfine in- 
teraction. The TDA, which also incorporates the same 
quasiparticle/constituent running masses as the RPA, in 
general provides a qualitatively comparable spin splitting 
description except for the ir-p difference. Of course by in- 
creasing the hyperfine strength it would be possible for 
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the TDA to account for the large ir-p splitting, however, 
this enhanced interaction would then generate an over- 
prediction for the iv'-p' and other splittings. For a com- 
prehensive description the TDA would most likely require 
a more complicated hyperfine interaction with tuning and 
in this sense shares the same difRculties as constituent, 
non-chiral models mentioned in the introduction. The 
attractive feature of the RPA-BCS approach is the abil- 
ity to obtain a good description with minimal parameters 
which is a common goal in all approaches to QCD and 
hadron structure. 

Application to the isoscalar hyperfine splitting is not 
possible without a proper 77-77' mixing calculation. How- 
ever, using the strange current mass of 25 MeV and the 
preferred potential, model 4, the splitting for the pure 
ss mesons is 300 MeV, corresponding to a pseudoscalar 
mass of 720 MeV and a meson at 1020 MeV. The ex- 
tracted strange constituent mass was 192 MeV. As the 
physical 77 and rj' masses are respectively 547 and 958 
MeV, mixing effects are significant and it will be of in- 
terest to see the importance of the hyperfine interaction 
in a rigorous mixing analysis. 

In addition to an improved meson spin spectrum, three 
of the model hyperfine interactions markedly increase the 
quark condensate which previous analyses 10, 16. 11^ IT^ 
predicted too low, around (qq) — —(110 MeV)'^. This 
was also noted in Ref. [22|. For models 2, 3 and 4 
the new condensate varies between —(165 MeV)"^ and 
— (200 MeV)'^ in the chiral limit, a noticeable improve- 
ment but still below accepted values spanning the inter- 
val between -(220 MeV)^ and -(260 MeV)^. The pion 
decay constant, also improves but only marginally. 
These shifts are in the correct direction, as first noted by 
Alkofer and Lagae [2^, but complete agreement in this 
model is not possible without generating very large self 
energies which distort the meson spectrum. Because the 
decay constant is a matrix element connecting the ground 
state (model vacuum), the low calculated values also re- 
flect shortcomings with the BCS vacuum. As detailed 
in Ref. 0| the use of the superior RPA vacuum signif- 
icantly increases /tt (although not quite to the physical 
value). It would therefore be interesting to repeat this 
hyperfine calculation with an improved vacuum. 

Naively, it is expected that the transverse potential 
should, as in the quantum mechanical quark model, de- 
crease the mass of pseudoscalar states. However, it is 
important to distinguish between level splitting and ab- 
solute level shifts. The hyperfine interaction does indeed 
provide a level splitting with difference proportional to 
hyperfine strength, but it also increases the quasiparti- 
cle self-energy and thus the effective constituent quark 
mass as well. Consequently, both pseudoscalar and vec- 
tor meson masses increase which then in turn requires 
a reduction in the current quark mass to reproduce the 
observed spectra. 

This reduction in current quark mass was also neces- 
sary to describe states in bottomonium. The reported 
but unconfirmed rjb state has a mass of 9300 MeV |27j . 



clearly below predictions (see Table HTjl which are closer to 
non-relativistic perturbative QCD (NRPQCD) and lat- 
tice calculations that predict a much smaller T-T^f, 
splitting of about 40 MeV with an error of about 20 
to 30 MeV. A recent paper ,2^J lowers this error to 10 
(th)lg {5a s) MeV reflecting uncertainties in theory and 
the strong coupling constant. The NRPQCD calcula- 
tions, in particular, suffer from uncertainty in nonper- 
turbative corrections that are usually parameterized in 
terms of condensates. Hence model calculations are still 
needed and the RPA predicts a splitting in rough agree- 
ment, but somewhat larger, than these theoretical expec- 
tations. Interestingly, the structure of the coupling (see 
Eqs. H31|l . H32|) and H33|) ) is such to increase the split- 
ting from the TDA value of around 20 MeV to the RPA 
prediction of 60 MeV. Although this is a minimal, sec- 
ondary effect when compared to the absolute mass scale 
involved, it is still important for the relatively small hy- 
perflne separation. This analysis therefore predicts that 
the 77f, meson mass should be around 9400 MeV and new 
results from spectroscopic studies at the B-factories are 
eagerly awaited. 



B. Comparison to other hadronic approaches 

It is instructive to make contact with alternative for- 
mulations and to further discuss the smaller current and 
constituent quark masses in our model. It is well known 
that a running current quark mass emerges in one loop 
QCD and values can be extracted at the perturbative 
Mz scale. Usually quoted, however, are their values at a 
much lower renormalization point such as 2 GeV in the 
MS scheme. While the "experimental" bare quark values 
listed in Table UTI are obtained from measurement, they 
rely on signiflcant model input either from chiral pertur- 
bation theory, as in the case of the u, d, s masses, or from 
heavy quark effective theory for the c, b quarks. They are 
thus not directly observable and entail both uncertainty 
and also ambiguity ^30]. It is therefore more appropriate 
to regard them as parameters in the Lagrangian (Hamil- 
tonian), subject to renormalization. Since our model 
kernel has been previously renormalized the only rem- 
nant of current quark mass running is the momentum 
dependence of the dressed quark mass Mq (p) . Effectively, 
the current quark renormalization point dependence has 
been converted into a constituent mass momentum de- 
pendence analogous to the Schwinger-Dyson treatment. 
Therefore when comparing to current masses from other 
approaches one should use the constituent running mass 
evaluated at the alternative model's current mass scale, 
e.g. Mq{p = 2 GeV), and not the smaller 777,^ parameter 
appearing in our Hamiltonian. Figure|21plots the running 
dependence of our model 4 dressed quark masses for dif- 
ferent flavors. As this flgure indicates at p = 2 GeV, the 
scaled effective current quark masses are still lower than 
in other approaches, however they are much larger than 
our bare Hamiltonian values and provide a more realis- 
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Running quark mass 

(model 4) 




contrasted with the estimate J = 0.48 using the physical 
TT, p, K and K* masses. The difference reflects the need 
to include dynamical fermions in the lattice calculations. 



p(GeV) 



FIG. 2: Constituent u/d, s, c and b running quark masses 
from Eq. fM . 



tic comparison. Related, the larger constituent masses 
in conventional quark models (ranges listed in Table II 
under "experiment") partially reflects missing dynamics 
from field theoretical self-energies that we explicitly in- 
clude. The validity of any quark approach should conse- 
quently be judged more by the robustness of observable 
prediction (e.g. spectrum) rather than specific quark 
values. This is the case in the present analysis as in- 
dicated in Table ^ where the resulting quark masses 
are small when compared to the typically quoted values 
but the predicted meson masses are reasonable. We sub- 
mit these smaller current quark values are representative 
of this simple, minimal parameter approach because the 
masses were uniformly small for all four, markedly differ- 
ent hyperfine interactions. An improved, rigorous treat- 
ment entailing a complicated combined quark-gluon sec- 
tor diagonalization of the exact hyperfine Hamiltonian, 
Eq. |(7J), is in progress which will firmly ascertain if small 
quark values are required. 

Finally, predictions for the "J parameter" 



J = M, 



dMv 
dMf^ 



(44) 



are presented which has been proposed as a reliabil- 
ity measure for quenched lattice computations of light 
hadron masses IsJ. Here Mps, My are the pseudoscalar, 
vector meson masses and Mn is the (reference) vector 
mass determined by the intersection of the line My ~ 
1.8 Mps with the plot of Mp^ versus My- If the vec- 
tor meson mass is linear in the current quark mass and, 
as indicated by Eq. ((2Jl, the pseudoscalar scales as the 
square root, then a sensitive lattice chiral extrapolation 
is not required to evaluate J and the attending errors can 
be avoided. Results from unextrapolated quenched lat- 
tice simulations [sj predicts J = 0.37 which should be 



RPA 




0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 
My (GeV) 

FIG. 3: Determining the J parameter. Tlie solid curve cor- 
responds to My ~ 1.8 AIps while the other lines are linear 
fits to tlie TDA and RPA predictions. 

Figure 13 shows the TDA and RPA masses and the 
curve My = 1.8 Mps- As anticipated, the RPA points 
scale linearly and extrapolate to zero pion mass at a 
vector mass of approximately 780 MeV. A linear fit to 
the RPA gives a reference mass of Mp — 880 MeV and 
J = 0.42, in more reasonable agreement with the esti- 
mate from data. Surprisingly, the TDA points also scale 
linearly even though they do not yield a zero mass pion 
in the chiral limit. This lowers the reference vector mass 
and produces a J parameter of 0.36. It will prove in- 
structive to confront these predictions with dynamical 
quark lattice simulations especially those using realistic 
sea quark masses. 



IV. CONCLUSION 

As a consequence of this study, the relative impor- 
tance of chiral symmetry and the hyperfine interaction 
is clearer: spin splittings in heavy quark systems are not 
governed by chiral symmetry and only require a hyperfine 
interaction. However, for light mesons chiral symmetry 
is important and is essential for describing the ir-p mass 
difference in a minimally parameterized, heavily con- 
strained model such as the one advocated here. Indeed 
the RPA-BCS many-body approach provides a reason- 
able description of the pseudoscalar-vector spectrum for 
both light and heavy mesons with a common Hamiltonian 
containing only the current quark masses and two dy- 
namical parameters. By explicitly incorporating this im- 
portant symmetry of QCD, a small pion mass is dynam- 
ically generated without the necessity of tuning a com- 
plicated hyperfine potential as typically done in conven- 
tional quark models. Furthermore, including a hyperfine 
interaction in this many-body approach improves both 
the pion decay constant and the quark condensate pre- 
dictions which previously have been calculated too low. 
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The hyperfine interaction also enhances the self energy 
contribution to the quark kinetic energy which necessi- 
tates using much smaller current quark masses. Lastly, 
the RPA J parameter is closer to data than quenched 
lattice results and it will be interesting to compare with 
dynamical quark lattice simulations. 

Future work includes reanalyzing the glueball, meson 
and hybrid spectra with the hyperfine potential and ex- 
amining other short range interactions such as the tensor, 
a^aj - a^aj, and higher dimensional terms from excluded 
Fock space components 23J . Extensions of this approach 
to highly excited hadron states will also be of interest, 
based upon the need js^ for new relativistic, chirally 
invariant models with a nontrivial vacuum. Finally, in- 
vestigations of baryons should also be fruitful as previous 



non- hyperfine calculations [33 only predict about half of 
the observed iV-A splitting. 
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